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Abstract. The luminosity distance - redshift relation is analytically given 
for generalized Randall-Sundrum type II brane-world models containing Weyl 
fluid either as dark radiation or as a radiation field from the brane. The 
derived expressions contain both elementary functions and elliptic integrals of 
the first and second kind. First we derive the relation for models with the 
Randall-Sundrum fine-tuning. Then we generalize the method for models with 
cosmological constant. The most interesting models contain small amounts of 
Weyl fluid, expected to be in good accordance with supernova data. The derived 
analytical results are suitable for testing brane-world models with Weyl fluid when 
future supernova data at higher redshifts will be available. 



1. Introduction 

At present the Universe is considered a general relativistic Friedmann space-time with 
flat spatial sections, containing more than 70% dark energy and at about 25% of dark 
matter. Dark energy could be simply a cosmological constant A, or quintessence or 
something entirely different. There is no widely accepted explanations for the nature 
of any of the dark matter or dark energy (even the existence of the cosmological 
constant remains unexplained). 

An alternative to introducing dark matter would be to modify the law of 
gravitation, like in MOND [T] and its relativistic generalization [2]. These theories 
are compatible with the Large scale structure of the Universe [3] . However in spite of 
the successes, certain problems were signaled on smaller scales [4]. 

Quite remarkably, supernova data, which in the traditional interpretation yield to 
the existence of dark energy, can be explained by certain f(R) [5] or inverse curvature 
gravity models [5]. However the parameter range, in which the latter is in goood 
agrement with the supernova data, also presents stability problems [7] . 

Modifications of the gravitational interaction could also occur by enriching the 
space-time with extra dimensions. Originally pioneered by Kaluza and Klein, such 
theories contained compact extra dimensions. The so-called brane-world models, 
motivated by string / M-theory, containing our observable 4-dimensional universe 
(the brane) as a hypersurface, were introduced in [8], [9] and [10], the latter model 
allowing for a non-compact extra dimension. 
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The curved generalizations of the model presented in [lOj have evolved into a 
5-dimensional alternative to general relativity, in which gravity has more degrees 
of freedom. In contrast with standard model fields, these evolve in the whole 5- 
dimensional bulk. In this generalized Randall-Sundrum type II (RS) theory, the brane 
has a tension A and gravitational dynamics is governed by the 5-dimensional Einstein 
equation. Its projections to our observable 4-dimensional universe (the brane) are 
the twice contracted Gauss equation, the Codazzi equation and an effective Einstein 
equation, the latter being obtained by employing the junction conditions across the 
brane [TT]. The effective Einstein equation (for the case of symmetric embedding 
and no other contribution to the bulk-energy-momentum than a bulk cosmological 
constant) was first given in a covariant form in [12j . Supplementing this by the pull- 
back to the brane of the bulk energy momentum tensor H a b, which is 

V ab = — {glg d b n cd ) (1) 

(with k 2 the bulk coupling constant and g a b the induced metric on the brane) the 
effective Einstein equation reads [XT] : 

Gab = -A.g ah + K 2 T ab + K A S ab - Sab + V ah ■ (2) 

Here k 2 is the brane coupling constant, related to the bulk coupling constant and the 
brane tension A as 6k 2 = k A, and 

A = y A - y n c n d U cd (3) 

represents a cosmological "constant" which possibly varies due to the normal 
projection of the bulk energy-momentum tensor (this includes the contribution —Ag ao 
due to the bulk cosmological constant A). The source term S a b is quadratic in the 
brane energy- momentum tensor T ab : 



i> ab - ~ 



-T ac T b c + l -TT ab ~ 9 -f (~T cd T cd + l -T 2 



(4) 



and £ ab is the electric part of the bulk Weyl tensor C a b c d, given as 

£ ac = C abc dn b n d . (5) 

In a cosmological context and suppressing any energy exchange between the brane 
and the bulk, this latter term generates the so-called dark radiation. Otherwise it can 
be called a Weyl fluid. 

A review of many aspects related to the theories described by the effective Einstein 
equation © can be found in [T3] . Both early cosmology [T3] and gravitational collapse 
[15]- [20] are essentially modified in these theories. There is also possible to replace 
dark matter with geometric effects in the interpretation of galactic rotation curves, 
weak lensing and galaxy cluster dynamics [2"T] . 

The possible modifications of gravitational dynamics are even more versatile in 
the so-called induced gravity models. These can be regarded as brane- world models 
enhanced with the first quantum-correction arising from the interaction of the brane 
matter with bulk gravity. The induced gravity correction couples to the 5-dimensional 
Einstein-Hilbert action with the coupling constant 7« 2 /k 2 . The simplest of such 
models, the DGP model was introduced in [25]. This model however suffers from linear 
instabilities (ghost modes in the perturbations), as shown for de Sitter branes [23]. The 
ghost modes withstand even the introduction of a second brane [24]. Generalizations 
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Figure 1. (Color online) A diagram presenting various brane- world models 
and their inter-relations. LWRS is the generalized Randall-Sundrum model with 
cosmological constant and a Weyl fluid reflecting a brane radiating into the bulk 
during nowadays or at least until recent cosmological times. 



of the DGP model are discussed covariantly in [25] and [26] when the embedding is 
symmetric, and in \27\ when it is asymmetric. In these models the role of the effective 
Einstein equation ([2]) is taken by a more complicated equation (see for example Eq. 
(29) of [27]), which contains the square of the Einstein tensor G a b- This implies that 
in certain sense the degree of nonlinearity of the theory is squared. In a cosmological 
setup the square root of this equation can be taken, leading to a set of modified 
Friedmann and Raychaudhuri equations, which however contain a sign ambiguity 
e = ±1 due to the involved square root. These are called the BRANE1 [DGP(-)] 
branch for e = -1 and BRANE2 [DGP(+)] for e = 1 in the terminology of [25] 
[or [28], respectively]. Both the original Randall-Sundrum type II model and the 
DGP model are contained as special subcases. Notably, the BRANE2 branch contains 
cosmological models which self-accelerate at late-times. We give in Fig Q] a diagram 
containing a classification of these theories and how they emerge as different limits 
from each other. 

In this paper we discuss analytically the luminosity distance - redshift relation 
in various generalized Randall-Sundrum type II brane-world models described by Eq. 
|2|). Our analytical approach can enhance the confrontation of these models with 
current and most notably with future supernova observations. We note that recently 
analytical results have been given in Ref. [29] for a wide class of phantom Friedmann 
cosmologies too, in terms of elementary and Weierstrass elliptic functions. 

In section 2 we review the notion of luminosity distance, its relation with the 
redshift and how these can be measured independently. This section was included 
mainly for didactical purposes. 

In section 3 we review the modification of this relation in the Randall-Sundrum 
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type II brane- world scenario. These include the introduction of the parameters fl\ and 
fid which can be traced back to the source terms S a b and £ a & of the modified Einstein 
equation @. The other cosmological parameters are fl p , representing (baryonic and 
dark) matter and f^A- We do not include bulk sources in the analysis, with the notable 
exception of a bulk cosmological constant. 

Section 4 contains the derivation of the analytic expression for the luminosity 
distance - redshift relation for the brane-worlds which are closest to the original 
Randall- Sundrum scenario [10], thus with no cosmological constant (Randall-Sundrum 
fine-tuning). The generic expression (|35[) of the luminosity distance derived here is 
given in terms of elementary functions and elliptic integrals of the first and second 
kind. From this most generic case we take the subsequent limits: fid = (subsection 
4.2), fix = (subsection 4.3); and both fid = fix = 0, this being the general relativistic 
Einstein-de Sitter case (subsection 4.4). 

Such models however could not allow for late-time acceleration, therefore in 
section 5 we discuss the luminosity distance - redshift relation for brane-worlds with A. 
First we present in subsection 5.1 a class of models, for which the luminosity distance 
can be given in terms of elementary functions alone. These models are characterized 
by an extremely low value of the brane tension, thus are in conflict with various 
constraints on brane-world models. 

Next, in subsection 5.2 we discuss brane-worlds for which the brane-characteristic 
contributions fix and fid represent small perturbations. This is a good assumption as 
observational evidences suggest that general relativity is a sufficiently accurate theory 
of the universe, and as such the deviations from it could not be very high, at least 
at late-times. We give analytical expressions in terms of both elementary functions 
and elliptic integrals of the first and second kind for the luminosity distance, to first 
order accuracy in the chosen small parameters of the model. Some of the most lengthy 
computations needed in order to achieve the result are presented in the Appendix. 

Section 6 contains the concluding remarks. 

Throughout the paper c = 1 was employed. 

2. The luminosity-redshift relation 

The Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 



describes a homogeneous and isotropic universe. Here r is cosmological time, (r, 9, if) 
are comoving coordinates, a is the scale factor and k = 0, ±1 the curvature index. 
The proper radial distance is defined as ar. A useful alternative form of the FLRW 
metric is 




dr 



,2 



dsi 



FLRW — 



dr 2 + a 2 (r) [d X 2 + r 2 ( X ] k) (dd 2 + sin 2 Otikp 2 )] 



(7) 



with 




k = 1 
k = 
k = - 



1 . 



(8) 



X being an other comoving radial coordinate 
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Figure 2. (Color online) A schematic representation of the propagation in the 
curved space-time of the light emitted by a supernova explosion in a distant galaxy 
and collected on the telescope mirror. A dimensional magnification (a la Wheeler) 
shows the elementary area dA normal to the direction of propagation n and the 
elementary solid angle dQ around n. 



If a photon stream emitted by an astrophysical light source travel without 
collisions, the number of photons d/V 7 from a comoving elementary volume of the 
6-dimcnsional phase space (x, p) is conserved |30j . Thus the phase space density 

1{,,p) ~d 3 £d s p Lo 2 dTdAdujdn {) 

of a photon stream is constant in time. Here u> denotes the frequency of the photons, 
dA and dfl stand for the elementary area normal to the direction of propagation and for 
the elementary solid angle around the direction of propagation, respectively (see Fig 
[2]). Eq. HJ holds true for any kind of cosmological evolution, provided d 3 x cx drdA 
and d 3 p oc uj 2 dujdtt are valid for the photons [30) . The luminosity of the source is 
C = dE em /dt em (total energy produced in unit time; the suffix em refers to emission). 

A telescope detects the photon flux T = dE rec / 'dT rec / 'Am (the suffix rec refers 
to reception). This is the energy detected during unit time on the telescope mirror 
surface Am- (The surface Am is understood to be perpendicular to the incident light 
stream.) 

From their definition, one can easily find a relation between T and C: 

J~ Am _ dE rec /dr rec (10) 
C dE em J ' dT ern 

As the energy of the photon stream in the comoving elementary phase space volume 
is dE = hu! dN, from Eq. © we find 



dE re r Am ( LOre.r\ dui rec dA rea dT r . 



dE ern At 



(11) 
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Here we have used that from the isotropy of the FLRW universe dQ rec — dQ em and 
we integrate the first to the solid angle encompassing the mirror surface, the second 
to the whole solid angle (cf. the definitions of E rec , E em ). In Eq. (1 1 1 1) A tot represents 
the proper area of a sphere centered in the light source and containing the reception 
point on its surface, at the time of reception. 

Due to cosmological evolution the elementary area dA changes as a 2 and the 
frequency of the light is redshifted during cosmic expansion, u a 1/a [30] . In the 
cosmological evolution of the comoving elementary phase space volume element dw 
changes accordingly: dm oc 1/a. Therefore 

L A tot \a-oJ 

where ao is the present value of the scale factor, and a is understood to be the scale 
factor at emission time. In the FLRW universe the proper area of a sphere with 
comoving radius r em is A to t = 47ra 2 l r 2 , T1 , an d the redshift z is defined as 

l + z=^L. (13) 
The luminosity distance oIl is defined as in Euclidean geometry: 

/ c \ 1/2 

dL ^'' = \Mr) = a O r em(l + Z) . (14) 

This definition is rigorous as long as we are dealing with the (homogeneous and 
isotropic) FLRW universe (irrespective of the value of the curvature index k) and 
the radius of a sphere is measured in the proper distance ra (the FLRW metric © 
guarantees that the surface of a sphere with radius ra is 47ra 2 r 2 ). 

According to Eq. |(7j) the comoving coordinate r em can be written in terms of an 
other radial comoving coordinate Xem (representing the location of the source): 

d L (z) = a (1 + z)H (xem] k) . (15) 

Disregarding possible deflections by perturbations of the FLRW universe, a light ray 
follows radial null geodesies of the FLRW metric, characterized by d\ = dr/a(r) = 
da/a 2 H. Here H — a/a is the Hubble parameter. Then 

f a " da f a ° da 

Xem=x{aem)= ~= 2 „ , T ■ (16) 

Ja Bm aa Ja am a H (a) 

By employing Eq. ()13|) the radial variable \ can a l so be expressed in terms of an 
integral over the redshift as 

/ \ 1 f" dz ' 

Xem{z) = — / "777 A ' ( 17 J 

a J H (z') 



which completes the definition (|15[1 of the luminosity distance 4 in terms of the 
redshift z. 

Differentiating Eq. (fT5")) with \ given by Eq. ([T7]) with respect to z gives 
1 



kd 2 L (z) ' 


- 1/2 d 




a 2 (l + z) 2 _ 


dz 


.1 + 2. 



1 ^LVJ /-|o\ 

H{z) L «o(l + ^) 2 J dz [l + z\ ' K ' 

therefore if independent measurements of dh and z are available for a set of light 
sources, the Hubble-parameter H (z) and in consequence the cosmological dynamics 
can be determined. 
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From the combined measurements of the large-scale structure of the Universe [31] , 
[32"] and of the structure of the cosmic microwave background [33J the conclusion was 
reached that the space geometry has flat spatial sections. Therefore in what follows 
we consider k = 0. Then the luminosity distance-redshift relation becomes 

d - {z) = {1+z) [-iw)- (19) 

In practice, the function d^(z) is conveniently measured with distant supernovae 
of type la. The luminosity is evaluated by photometry, while the redshift from 
spectroscopic analysis of the host galaxy. 

Each cosmological model has its own prediction for the shape of the function 
dL{z) [see Eq. (TT5]) with \ given by Eq. (fT?]) for generic fc, or Eq. ([1^]) for k = 0]. 
This is how the measured cLl(z) data turn into a cosmological test. 



3. The luminosity-redshift relation in Randall-Sundrum type II 
brane- worlds 

We consider FLRW branes with k = and brane cosmological constant A, embedded 
symmetrically. The bulk is the Vaidya-anti de Sitter space-time with cosmological 
constant A, and it contains bulk black holes with masses m on both sides of the 
brane. The black hole masses can change if the brane radiates into the bulk. An 
ansatz comparable with structure formation has been advanced in [34] for the Weyl 
fluid m/a 4 for the case when the brane radiates, m — moa a , where mo is a constant 
and a — 2, 3. For a — the Weyl fluid is known as dark radiation and then the 
bulk space-time becomes Schwarzschild-anti de Sitter. The brane tension and the two 
cosmological constants are inter-related as 

2A = k 2 A + k 2 A. (20) 

The Friedmann equation gives the Hubble parameter to A, m, the scale factor a 
and the matter energy density p on the brane: 

2 A n 2 p ( ^ 9 \ , 2m 



« =3 + tH 1 + ^ + ^- < 21 » 

In the matter dominated era the brane is dominated by dust, obeying the continuity 
equation 

p + 3Hp = 0, (22) 
which gives p ~ a~ 3 . We introduce the following dimensionless quantities: 

&tot = n A + n p + n x + n d , (23) 
n p = — ^ . n x = — ^ , (24) 

n d =_^° fi A = A_- (25) 

The subscript denotes the present value of the respective quantities. In terms of 
these notations the Friedmann equation becomes 



= Hi 



(26) 



The luminosity-redshift relation in brane- worlds: 



I. Analytical results 8 



In particular at present time this gives f2 tot = 1. Then the radial coordinate (|16p 
becomes 

Xem = -ft / -rjK ■ (27) 

J aem [n A a e + fl p a 3 a 3 + n d a^ a a a + 2 + n x a%] 1 

This is a complicated integral, which cannot be computed analytically in the majority 
of cases. In what follows we will analyze various specific cases of the above integral, 
when an analytic solution is possible. The cases a — 2, 3 represent the Weyl fluid 
compatible with structure formation, while a = represents the dark radiation. 

4. Branes with Randall- Sundrum fine-tuning 

In the original Randall-Sundrum scenario the bulk cosmological constant A is fine- 
tuned with the brane tension A such that cf. Eq. (|20|) the brane cosmological constant 
vanishes. For simplicity we also assume throughout this section a — 0. By imposing 
a vanishing cosmological constant on the brane, Q\ — such that the polynomial of 
rank 6 in the denominator of the integrand in Eq. (I27|) shrinks to a polynomial of rank 
3. Therefore its roots can be found analytically. Following general procedures, the 
luminosity distance - redshift relation can be then given analytically in terms of elliptic 
functions. This is done in the following subsection. In the second and third subsections 
of this chapter we discuss the limits — > (when the bulk is anti de Sitter) and the 
late-time universe limit p/X — > 0. The general relativistic (Einstein-deSitter) limit is 
found in the fourth subsection, when further fix — > is taken. 

4-1- Schwarzschild-AdS bulk 

With no brane cosmological constant, Eq. (f2"7| becomes: 



Xt 



a iiQiip 



ada 



n>« 2 + Sr«o 



1/2' 



(28) 



Following the method given in [35] we find the following roots of the denominator: 
a = — — 1 + 2 cosh — 



3n p \ 3 



= ^ (- 1 + cosh - + i V3 sinh - ) , (29) 



3n p V 3 3 

and its complex conjugate 0*. The auxiliary quantity ^ is defined as 

cosh* = i + __^. (30) 
We introduce the following real combinations of the complex roots 
6 1 = ^ = ^f-l + cosh^ 



2 3ft p V 3 



ai = ; — = —j= — Sinn — . (31) 



% V3ft p 3 
Then Eq. (|28|) is written conveniently as 

1 f a " ada 



Xern — 



3/2 jj l/2 , - „ 

H oUp J «™ ( a -a) (a-hf + al 



1/2' 



(32) 
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The integration can be carried out by employing the formulae (239.07) and (341.53) 
of Ref. m. We obtain 



1 



3/2 „ l/2 
a Jrt il p 



g(a + A) [F (</> , e) - F (ip em , e)] 



2gA [E (<f ,e)~E((fi em ,e)} 
sin ipo \/\ — e 2 sin 2 <^q 



sin (^emVT 



1 + COS ifio 



e 2 sin 2 i^em 



1 + COS tp e 



(33) 



where F (</?, e) is the elliptic integral of the first kind; E (ip, e) is the elliptic integral 
of the second kind (with variable ip and argument e); and we have introduced the 
following standard notations, cf. Ref. [55] : 

A 2 = (bi-af + a 2 , 
2 A + bi — a 



2A 



1 



.4 



A + a — a 



(34) 



A — a + a 

By employing Eqs. (fT3|) . (fl~7|) and (|19|) . after a lengthy, but straightforward calculation 
the luminosity distance-redshift relation emerges: 



Xd _(l + z)Q d 



1/2 



{B 1 +B 2 )[F (<p ,e)-F(<p em ,e)] 



2B, 



E{ip em ,e) -E(cp ,s) 



sin ipo 



vT 



e 2 sin 2 tpo sin tp em 



S 2 Sin 2 (fi em 



with 



Bi 



1 + COS ifio 

-3 1 / 4 (1 + 2 cosh I) 



1 + COS (f e 



(35) 



3 (4 cosh 



2 « _ ^1/4 ' 



£? 2 = 



4 cosh 2 f - 1 



1/4 



(36) 



and 



1 V3 cosh 
— I ^ 

2 2 A /4cosh 2 f -1 



(l + «) «d 



arccos < 



(i + z) n d 



12cosh 2 f-3-l-2coshf 



12cosh 2 f-3 + l + 2coshf 



-3fL 



..(37) 
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Here tp runs in the range 0..7r/2. In computing for other values of tp, we can use the 
following addition rules for the elliptic integrals: 

F(mn±tp,s) = 2mK±F(tp,s) , 

E(m%±tp,s) = 2mE ± E (tp, e) , (38) 

where K and E are the complete elliptic integrals of the first and second kind. Eqs. 
(|30| and ([35|) -([37 ]) represent the analytical expression of the luminosity distance- 
redshift relation for FLRW branes with Randall-Sundrum fine-tuning. They are 
given in terms of the well-known elliptic integrals of first and second kind, and the 
cosmological parameters ft p , Q\ and ft^. 



4-2. Limit of no black hole in the bulk 



In this subsection we consider the case ft^ = 0. The derivation follows closely the steps 
of the previous subsection, however the formulae are simpler. The auxiliary expression 
(|30[) for \P is well defined only for ft^ ^ and we have to address the question how to 
obtain suitable limits of the results derived for ft<; ^ 0. For any ft^ -C 

27ft A ft 2 



cosh 'I' 



But 



cosh ^ — cosh — 
3 



» 1 



4 cosh 2 — - 3 



4 coslr — 



as cosh(*/3) > 1 also holds. Thus 



cosh ■ 



± sinh • 



l/3 n 2/3 



3ft 1 / ^ 



(39) 



(40) 



(41) 



3 3 2ft d 

By employing Eq. (|4ip in the generic expressions derived in the preceding subsection, 
we obtain the luminosity distance-redshift relation in a very similar form to Eq. (|35p . 
but with different coefficients: 



2^3 (l + z) ft 1/6 




[F(ip ,e)-F(ip em ,s)] 



sin ip e 



e 2 sin 2 >p em 



where 



V3 
4 



<P 



COS <po 



l)ftf (1 



COS if e 



(42) 



ft 



1/3 



(^+l)ft^ /3 (l 



ft. 



1/3 



(43) 



Again, ip for this case emerges in the limit fid — > from the generic expression Eq. 
(|37[) . by employing Eq. ( |4"T|) as in the limiting process expressions of the type oo x 
appear. 
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4-3. Late-time universe limit 

In the late-time universe p <C A and in consequence fl\ = can be safely assumed. 
We keep however the dark radiation in the model. Eq. ([77]) simplifies considerably, 
and a straightforward integration gives the luminosity distance - redshift relation 



^ {z) = ^§7 (y/(n p + n d )(i + z) yjn p + n d (i + z)^ . (44) 

We can also prove that this result emerges as the il\ — > limit from the generic 
results, Eqs. {3D]) and ((35 ]) -(|37 ]> . When n x -> Eq. (gQ) gives * -> 0. Then Eqs. 
([57]) and ([II give 

e 2 =1 

V = arccos \ - 7 — — g — I , (45) 



2(i + ^)o d + n p 

B 2 = 1 = -Si . (46) 
By noting that -E (</?, 1) = siny, we obtain from Eq. (|35[) : 



. _ 2(l + z)»f 



0"p 



sin (^o sin ip e 



(47) 

1 + COS IfQ 1 + COS </? em 

By inserting the values ip em — tp (z) and po — tp (0), we recover the luminosity distance 
- redshift relation (|4"4"| . 

4-4- General relativistic (Einstein-de Sitter) limit 

The general relativistic limit of the luminosity distance - redshift relation for dust 
matter and k = = A (Einstein-de Sitter model) can be obtained by direct integration 
ofEq. (J27|: 



dl R (z) = 2 -£^[VT+-z-l} , (48) 

It is straightforward to check that the above result stems out from Eq. (|44|) by simply 
switching off the dark radiation. 

The general relativistic limit of the luminosity distance - redshift relation should 
also emerge in the limit Vl\ — ► of Eq. (|4^|) . To see this, we note that when £l\ — > 0, 
both p — > 7r and ipo — > 7r. Therefore the elliptic integrals of the first and second kind 
both tend to finite values, thus the differences evaluated at p> and ipo vanish. Then 
the only terms which should be carefully investigated are the last two terms of Eq. 
(|42p . which are of the type 0/0. By employing Eq. (|4"5)) . for the last term we obtain: 

r i/6 sin <PemV 1 - £ 2 sin 2 cp em ttl /6 

um »V ; = ., ,, , = ■ (49) 

Oa-o a 1 + cosyw 3 1 / i VTTl 

Accordingly, the second to last term gives 

lim ni/esinWl-^ai^go = 

fi^o A 1 + cos^o 3 x /4 v ; 

Adding everything together, we recover the general relativistic result (|4*5|) . 
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5. Branes with A 

In this section we discuss certain cases of Randall-Sundrum type brane-worlds with 
cosmological constant, for which analytical expressions for the luminosity-redshift 
relation can be found. 

5.1. A brane with analytically integrable luminosity distance-redshift relation 

If we do not impose the Randall-Sundrum fine-tuning in Eq. (|20[) and we keep the 
brane cosmological constant A, the polynomial in the denominator of the integrand in 
Eq. (|27[) can be simplified for certain values of the dimensionless Sl-s. In particular, 
if we choose 

n d = and 4:fl x n A = Q 2 p , (51) 

the expression under the square root of denominator becomes a quadratic expression, 
and the integral can be given in terms of elementary functions [37 : 



2 1 / 3 (l + z) f (l-h + h 2 ) [l + fc (! + «)] 



2 



4= K A/2 = - r,: ~l k in 



6# ttp /3 ^A 6 1 (1 + hf [l - h (1 + z) - h 2 (1 + zf 

r ,- 2(l + z)-h) , . 

+ 2 V3 arctan —= 2 V3 arctan — ■=!■ ) , (52) 

V3h V3h J 

with ft = (f7 p /2ft A ) 1/3 - 

The first condition (|5ip merely simplifies the bulk to an anti de Sitter space-time. 
The second condition (|51[) by contrast, yields to a much more serious constraint: 

A = 2A (53) 

The second condition (|51| , together with the constraint ([23l ) leads to a quadratic 
equation for For fi p = 0.27 this has two solutions [37] : 

A = 0.704, ft A = 0.02 6 (54) 

corresponding to the brane tension^] Ai = 38.375 x 10~ 60 TeV 4 and 

n A = 0.026, n x = 0.704 . (55) 

corresponding to the brane tension A2 = 1.4173 x 10 _60 TeV 4 . 

It is interesting to note that while solution (|55[) is ruled out by the recent 
supernova data, solution (|54"| is quite close to the present observational value of f^A 
|38j . From a brane point of view, however the value of the brane tension in the model 
([54]) is far too small, thus it does not describe our physical world. Indeed, all lower 
limits set for A are much higher than A2. 

In the two-brane model of Ref. [3] the minimal brane tension depends on the 
value of the Planck mass Mp and on the characteristic curvature scale of the bulk I 
as A m in = SMp/Aitl 2 9 . Table-top experiments on possible deviations from 
Newton's law currently probe gravity at sub-millimeter scales. As a result they 
constrain the characteristic curvature scale of the bulk to I < 44 /im. The brane 
tension therefore (in units c = 1 = K) is constrained as A > 715.887 TeV 4 . (For a 
detailed discussion see section 6 of [ID], where a slightly lower bound for the brane 
tension was derived, based on the previously available estimate I < 0.1 mm for the 

\ All values of the brane tension given in this subsection are in units c = 1 = h. 
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characteristic curvature scale of the bulk.) Big Bang Nucleosynthesis constraints give 
a much milder lower limit, A > 1 MeV 4 [H] . An astrophysical limit A > 5 x 10 s 
MeV 4 (depending on the equation of state of a neutron star) has also been derived 
|42j . This latter value of A m i n is in between the two previous lower limits. 

The interpretation of the model (|54[1 is the following. The condition (|53| on the 
models with small brane tension implies A = 0, thus the bulk becomes flat. As such, it 
has no effect on dynamics and the fifth dimension becomes superfluous. In fact what 
we face here is a GR model with stiff fluid scaling as a~ 6 . 

5.2. Branes with <C 1 and Q\ <C 1 

In this subsection we assume that both fl\ and fi^ are small, however we allow for 
arbitrary values of These assumptions are motivated by observational evidence 
that at present our universe is extremely close to a ACDM model. A Taylor series 
expansion of Eq. (|27p gives, to leading order in the small parameters: 

d AX d = ^ACDM + ^ a/a + ^ i (56) 



with 



d 



ACDM 



a Q (l + z) [ a ° da 



a 7 (l + z) f a ° da 



2H o Ja cm a 7 / 2 [n A a* + n p a 3 } 3/2 ' 
(Q) _ al a (l + z) [ ao a^da 

d " 2^„ J aem [nAa 3 + npa 3f 2 ■ [ > 

The first expression is the general relativistic luminosity distance - redshift relation in 
the presence of a cosmological constant (in the ACDM model) . The next two integrals 
represent the correction functions scaling the small coefficients ft\ and f^- 

All integrands have the same expression fl^a 3 + fl p aQ in the denominator. The 
roots of this cubic polynomial are: 

/<v 1/3 

^ a i(kT\ i+iVE ) (58) 

and /3*. Then d^ CDM can be rewritten as: 

,acdm = a o (1 + z ) r° da 

L rrnl/2 / 1/2 / sl/2/ ON 1/2/ n*^/ 2 ' ^ ' 

iio"A a a < {a — a) 1 (a — p) ' (a — p*J ' 

The integration can be carried out by employing Eq. (260.00) of [36j and we obtain 
the result: 

.acdm {z) = (l + z)[Fto,e) -F(<p,e)] 

with the variable <p and argument e of the elliptic integral of the first kind F (tp, e) 
given by 

2 1 , V3 
£ =2 + -' 
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arccos ■ 



i + *) ny 3 



i + Vs) n 1 / 3 



1/3 
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(61) 



(Note that e 2 is the same as in the case f?A = = Qd, while ip is different. Here 
< (f < 7r/2 while for other values of 93, we use Eqs. (I3"8l) .) 

It is relatively easy to integrate the contribution of the term linear in £l\ in terms 
of the variable t = a 3 / 4 . After a partial integration meant to reduce the powers in the 
denominator we employ 

,1/3 

[F(<po,e)-F(<p,e)] 



and obtain 



da 



ni 



,1/2 



1/2 



3 1 /4 ao ^ /3 



(62) 



15H W D 



ft, 



,1/2 



+ 



8fi A + 3f7 p (l + z)- 



fiA + fi p (l + z)- : 



1/2 



8r>^ /6 (1 + z ) 



[F(cp ,e)-F(cp,e)} , 



(63) 



with the variable <p and argument e given in Eq. (|61l) . 

The last term of Eq. ([56"]) is much more complicated to evaluate. For a = 1 
and 4 the source term ild merely contribute to ft p and respectively. The more 
interesting cases are for a = , 2, 3. The last term of Eq. (55) for a = 2 consits of 
elementary functions: 

1 + Z \ +Z , (64) 

3 3QpWflA + ftp 



(2) 



while and 1^ are more complicated to evalute, and we give details of the 
derivation in the Appendix. By passing to the variable z instead of a, we obtain: 



r (0) 



(1 + *) 



Afl A + 3ft c 



4^ A + 3(l + z) 3 fip 



+ Q p (i + z)y / r> A + (i + z) 3 «p 



,(l + z)Q 

3H n 5 p /3 



1/6 



-Id (<p, e) 



(65) 



and 



r(3) 



v /n ^ (i + z) v / ^ A + (i + z) 3 r!p 



2(l + z) 

3H Q n A /e n 2 p /3 



(66) 



where 



V3) 



:{F(ip ,e)-F(ip,e) 



^3(1 + V3)' 
- (3 + V3) [E(<po,e)-E(<p,e)] 



Vs) 
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x 



sin ipo\/ 1 — e 2 sin 2 ipo sin ipy/l — e 2 sin 2 tp 



} • (67) 



1 + (2 + V3) cos tpo 1 + (2 + V3) COS if 



Thus, the analytic expression of the generic luminosity distance - redshift relation on 
branes with cosmological constant and small values of fix and fid is given to first order 
accuracy in these small parameters by Eqs. (fS"6")) . (fBT))) , (j()3"|) - (f6"T)) . 

6. Concluding remarks 

The main purpose of this paper was to present the analytical formulation of the 
luminosity distance - redshift relation in the generalized Randall- Sundrum type II 
brane-world models containing a Weyl fluid either in the form of dark radiation or 
as radiation leaving the brane and feeding the bulk black holes. We have given the 
luminosity distance in terms of elementary functions and elliptical integrals of first and 
second type and we have also shown how the different limits arise from the generic 
result. Our results hold for: 

(a) Models with Randall- Sundrum fine-tuning (A = 0), with or without dark 
radiation from the bulk and with or without considerable contribution from the energy- 
momentum squared source terms, discussed in section 4. 

(b) The models discussed in subsection 5.1, obeying A = k 2 X/2, integrable in 
terms of elementary functions and 

(c) Models with a brane cosmological constant, discussed to first order accuracy 
in both the Weyl fluid and energy-momentum squared sources. 

This last class of models, presented in subsection 5.2 in the latest times of the 
cosmological evolution are only slightly different from the A CDM model, as they have 
fid <C 1 and fl\ <C 1. The derived modifications in the luminosity distance - redshift 
formula then represent corrections to the corresponding formula of the ACDM model. 

While the focus of the present paper is the integrability of the luminosity distance 
- redshift relation in various brane-world models, in a forthcoming paper [43] we 
will discuss how well the presently available supernova data support the brane-world 
models with a small amount of Weyl fluid. 
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Appendix A. The evaluation of the integral Id 

We can integrate the last term of Eq. (55) for a = 0, 3, as follows. First we pass to the 
variable t = a?l 2 and we perform a partial integration in order to reduce the powers 
in the denominator of the integrand 
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and 



(i + ^i /2 

3H a a fl 2 



At 2 + 3%^ 



■(i+£) »j / 2 r° 



^2 _|_ ^p°o 



r(3) _ 



3/2 



f 2 



3/2 



(A.l) 



(1 + *) 



2(l + z) 
9i/oao« A /2 "p 



t 2 / 3 dt 



£2 _|_ ^P a 



(A.2) 



By a change of the integration variable to x = Q^^t 3 ^ /Q^J^aa we can employ Eq. 
(260.52) of [35] in order to evaluate the remaining integral: 



t 2 / 3 dt 



3/2 



j-2 _|_ ^p a o 



n A 



2n 



1/3 



1/3 
A 
Q l/3 

P 
,V3 
\ 



,1/3 
V "0 



x 2 iix 



1/3 



1/3 



(A.3) 



Here (y>, e) is given by Eqs. (167 
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